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ABSTRACT 

We analyze extensive data on NGC 1377 from the Spitzer Infrared Nearby 
Galaxies Survey (SINGS). Within the category of nascent starbursts, that we 
previously selected as having infrared to radio continuum ratios in large excess of 
the average, and containing hot dust, NGC 1377 has the largest infrared excess 
yet measured. Optical imaging reveals a morphological distortion suggestive of 
a recent accretion event. Infrared spectroscopy reveals a compact and opaque 
source dominated by a hot, self-absorbed continuum (r ~ 20 in the lOyum sili- 
cate band). We provide physical evidence against non-stellar activity being the 
heating source. HII regions are detected through the single [Nell] line, probing 
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< 1% of the ionizing radiation. Not only is the optical depth in different gas 
and dust phases very high, but > 85% of ionizing photons are suppressed by 
dust. The only other detected emission features are molecular hydrogen lines, 
arguably excited mainly by shocks, besides photodissociation regions, and weak 
aromatic bands. The new observations support our interpretation in terms of an 
extremely young starburst (< IMyr). More generally, galaxies deficient in radio 
synchrotron are likely observed within a few Myr of the onset of a starburst and 
after a long quiescence, prior to the replenishment of the interstellar medium 
with cosmic rays. The similar infrared-radio properties of NGC 1377 and some 
infrared-luminous galaxies suggest that NGC 1377 constitutes an archetype to 
better understand starburst evolution. Although rare locally because observed 
in a brief evolutionary stage, nascent starbursts may represent a non-negligible 
fraction of merger- induced starbursts dominating deep infrared counts. Since 
they differ dramatically from usual starburst templates, they bear important 
consequences for the interpretation of deep surveys. 

Subject headings: dust, extinction — galaxies: individual (NGC 1377) — galax- 
ies: ISM — galaxies: starburst — galaxies: evolution — infrared: galaxies 



1. Introduction 

The infrared-radio correlation of star-forming galaxies (Helou et al. 1985; de Jong et 
al. 1985) suffers very few exceptions. At centimeter wavelengths, the radio continuum of 
massive galaxies is dominated by synchrotron emission, from cosmic ray electrons previously 
accelerated in type-II supernova remnants, propagating in the interstellar magnetic field, 
and decaying in less than 10^ years (Condon 1992). The emission from thermal electrons 
in HII regions usually contributes only of the order of 20% of the total emission of galaxies 
at 3 cm, and of the order of 5% at 20 cm. The far-infrared continuum measures the energy 
peak of the dust emission, and is a fair tracer of the instantaneous star formation rate for 
dusty starburst galaxies. In more quiescent galaxies, a significant fraction of the far-infrared 
emission can originate from dust heated by intermediate-mass stars ( Lonsdale- Persson & 
Helou 1987; Sauvage & Thuan 1992). 

NGC 1377 is the most extreme example known so far of the rare category of non-dwarf 
galaxies which have infrared-to-radio ratios in significant excess of the average, and where 
the infrared excess is not attributable to cirrus-like emission ; this is because the high dust 
temperatures of the selected galaxies imply high radiation field intensities (Roussel et al. 
2003, 2005). In NGC 1377, synchrotron emission is deficient by at least a factor 37 (i.e. by 
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> 8a) with respect to normal galaxies. In addition, HII regions were not detected through 
near-infrared hydrogen recombination lines nor through the thermal radio continuum, and 
the upper limits that we obtained imply that they are fainter than expected for normal 
starbursts by at least 70%, likely due to the fact that most ionizing photons are absorbed 
by dust. 

The infrared and radio emission both have their source in the star formation activity. 
Given the lack of direct connection between the physical processes regulating both types of 
emission, however, the infrared-radio correlation can hold only if strong coupling mechanisms 
operate (Beck & GoUa 1988; Helou & Bicay 1993; Niklas & Beck 1997). First, the production 
rate of cosmic rays has to be roughly proportional to that of dust-heating photons (which is 
achieved for a constant initial mass function). In addition, the magnetic field intensity has 
to be coupled with either the gas density or the star formation rate (e.g. turbulence driven 
by massive star formation inducing rapid amplification of the magnetic field). Empirically, 
starbursts seem to constantly adjust their magnetic field so that the ratio of magnetic energy 
density to radiation energy density remains constant, including ultraluminous galaxies with 
very intense radiation fields (Condon et al. 1991; Lisenfeld et al. 1996). 

The infrared emission and the radio continuum both persist for long periods subsequent 
to an instantaneous episode of star formation, but follow very different timescales. In partic- 
ular, cosmic rays are accelerated only about 4 Myr after the birth of massive stars evolving 
into type-II supernovae, whereas dust is heated almost instantaneously. This time delay 
implies that the relative amounts of infrared and radio emission are expected to vary signif- 
icantly, which can in principle be exploited to constrain the age of individual star forming 
regions, i.e. regions with relatively simple star formation histories. Galaxies, however, host 
young stellar populations which span an extended age range, and star formation generally 
goes on for long periods, either continuously or in the form of recurrent bursts. For large 
regions within galaxies, systematic variations of the infrared-to-radio ratio are thus in prac- 
tice erased by complex star formation histories. Non-uniform interstellar medium properties 
(structure and dust content) will also induce further variations decoupled from age. 

Discrete sources with radio spectra dominated by optically-thick thermal emission, from 
ultracompact HII regions, have been found and studied at high angular resolution within 
starbursting galaxies, most of them dwarfs: e.g. by Turner et al. (1998) in NGC5253, 
Kobulnicky & Johnson (1999) in He 2-10, Tarchi et al. (2000) in NGC2146, Beck et al. 
(2000) in several Wolf-Rayet galaxies. Beck et al. (2002) in IIZw40. These systems thus 
contain a number of very young star forming sites. On global scales, the radio continuum of 
blue compact dwarf galaxies can be dominated by the thermal component. In such objects, 
both the radio and the infrared emission deviate significantly from the scaling relations with 
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star formation rate followed by more massive galaxies (Hunt et al. 2005, and references 
therein) : the radio continuum is depleted because of the delayed injection of cosmic rays in 
a young starburst, and possibly because of enhanced cosmic ray escape ; and the infrared 
emission is depleted because of lower dust content, different dust properties, and higher 
porosity of the interstellar medium. 

We selected the galaxies that we term nascent starbursts from an infrared flux-limited 
sample by their high infrared to radio flux ratios and their high dust temperatures (Roussel 
et al. 2005). These criteria are designed to select galaxies hosting very young starbursts, like 
the galaxies mentioned above, but constituting a very different population: their infrared 
luminosities are much higher, they are more massive, and if cosmic ray escape is less efficient 
than in dwarfs, the radio deficiency also implies that cosmic rays produced by past star 
formation episodes have decayed, and thus have been injected more than ~ 100 Myr ago. In 
other words, the intensity contrast between the current starburst and the past star formation 
integrated over several tens of million years is very high. While classical starburst galaxies 
also contain very young star forming sites, our view of these is confused by recurrent burst 
episodes or continuous star formation in the recent past, which produce average signatures 
typical of evolved bursts, mimicking single populations of ages ~ 5-10 Myr (e.g. Calzetti 
1997). Nascent starbursts, on the other hand, offer an ideally simplified view because of 
the dominance of a single and quasi-instantaneous burst of star formation (previous stellar 
generations being too old to contribute to the production of star formation tracers). 

The nascent starbursts are extremely rare. They represent a negligible fraction {^ 1%) 
of an infrared fiux-limited sample in the local universe such as the IRAS Faint Galaxy 
Sample, and about 16% of galactic systems whose flux density is higher at 60 /xm than at 
100 /xm (Roussel et al. 2003). These objects nevertheless deserve careful study, because they 
constitute ideal settings to better understand the initial conditions and early evolution of 
starbursts, as well as the regulation of the infrared-radio correlation. 

In such objects, where the active regions are very young and compact, the most useful 
diagnostics on the nature of the heating source are expected to arise from the dust and the 
molecular gas phases, because ultraviolet radiation will be confined to dense ultracompact 
HII regions (Habing & Israel 1979). The mid-infrared spectral range is of particular interest 
not only to constrain the nature and excitation of the small-size dust grains and aromatic 
compounds, but also because it contains a series of forbidden and molecular hydrogen lines 
providing unique diagnostics on HII regions, photodissociation regions, and the outer layers 
of molecular clouds, and which are much less affected by extinction than the other diagnostics 
accessible in the optical and near-infrared domains. An ISOPHOT mid-infrared spectrum of 
NGC 1377 between 2.5 and 4.8 /xm and between 6 and 11.5 /xm was presented by Laureijs et 
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al. (2000), as part of a small sample of galaxies whose infrared flux density distribution peaks 
at 60;um. This spectrum appeared very peculiar in comparison with normal galaxies, con- 
taining a broad emission feature between 6 and 8.5 /xm in place of aromatic bands universally 
found in metal-rich star-forming galaxies, but lacked sufficient wavelength coverage to draw 
conclusions regarding the nature of the source and the excitation mechanisms of the dust 
and gas. NGC 1377 was selected to be part of the Spitzer Infrared Nearby Galaxies Survey 
(SINGS) sample (Kennicutt et al. 2003) to extend the range of explored -Fgo^mZ-^iooMm ratios, 
indicative of average dust temperatures. In this paper, we present the SINGS broadband 
imaging data (3 to 160 yum) and spectral imaging data (with continuous coverage between 
5 and 38;um) of NGC 1377, as well as ancillary optical imaging and spectroscopic data ob- 
tained for SINGS. We analyze them in view of constraining the geometry of the gas and 
dust reservoir, the excitation of the different phases of the interstellar medium probed in the 
5-38 /xm range, and what these properties imply for the nature of the activity in NGC 1377. 
The hypothesis of non-stellar activity is discussed in Section 9. 



2. SINGS data 

2.1. Spitzer spectral maps 

Low spectral resolution cubes (A/AA ~ 60-130) between 5 and 38/im and high spectral 
resolution cubes (A/AA ~ 600) between 10 and 37 /im centered on the nucleus of NGC 1377 
were acquired with the IRS instrument (Houck et al. 2004a) in February 2005 and August 
2004, respectively. For each spectral resolution mode, the full wavelength range is covered in 
two parts by two different slits: at low resolution, SL (5-14.5 /im) and LL (14-38 /im); and, 
at high resolution, SH (10-19. 5/im) and LH (19-37;um). Quadruple coverage is obtained at 
each point of the maps, except at the edges, by moving the slit in increments of half the slit 
size in each direction (see Kennicutt et al. (2003) and Smith et al. (2004) for details on the 
observing strategy in SINGS). Given that the infrared emission of NGC 1377 is not spatially 
resolved (Section 3), all the emission is recovered within the spectral maps, the smallest of 
which, in the SH module, is 23" x 15". The data were pre-processed with the S12 version of the 
Spitzer Science Center pipeline, and then assembled in spectral cubes with Cubism. Cubism 
is the software developed to rectify, project and co-add SINGS data, with the capability 
to subtract a spectral cube of the background and foreground (zodiacal) emission extracted 
from the data themselves, to mask pixels exhibiting an abnormal responsivity, and to perform 
spatial and spectral extractions (Smith et al., 2006, in preparation). In order to improve 
the signal to noise ratio of the high-resolution spectral data, an extra step was included in 
the processing of these data, outside of Cubism. In the SH and LH modules, the slit is 
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very narrow, only 6 pixels-wide on average, and the pixels with an abnormal responsivity 
are more numerous than in the SL and LL modules. Because of the minimal redundancy, 
masking pixels can create spurious features in the final cube. To mitigate this effect, the 
stacked diffraction images were inspected for obvious deviant pixels (i.e. features that are 
narrower than the spectral and spatial resolution), and these pixels were interpolated in each 
pointing along the wavelength dimension. This treatment is appropriate for infrared point 
sources like NGC 1377, with abrupt spatial variations and a smoothly-varying continuum, 
and leaves intact real spectral features. 



2.2. Spitzer broadband imaging 

Images in the four bands of the IRAC camera (Fazio et al. 2004), centered at effective 
wavelengths of 3.6, 4.5, 5.7 and 7.9 /xm (usually designated by their fiducial wavelengths 3.6, 
4.5, 5.8 and 8.0 /im), were obtained in August 2004. Scan maps in the three bands of the 
MIPS instrument (Rieke et al. 2004) at effective wavelengths of 24, 71 and 156 /xm (fiducial 
wavelengths of 24, 70 and 160 /im) were acquired in September 2004. The observing strategy 
and data reduction are described by Kennicutt et al. (2003). The full width at half maximum 
of the point spread function (PSF) is close to 2.0" at 7.9 /im, or 200 pc at the distance of 
NGC 1377 (21Mpc). Flux calibration uncertainties are of the order of 10% in the IRAC 
bands, and 5%, 10% and 15% in the MIPS 24, 71 and 156 /im bands. 

Broadband infrared flux densities of NGC 1377 are listed in Table 1. In the IRAC 
short-wavelength bands, both the total fluxes and the nuclear fluxes are provided. In all 
IRAC bands, the nuclear fluxes were measured within the smallest point-source aperture of 
diameter 4.88", with a background annulus of internal and external diameters of 4.88" and 
14.64", in which case aperture correction factors ranging from 1.213 to 1.584 are necessary 
(IRAC handbook). Use of the other point-source apertures tabulated in the IRAC handbook, 
for the 5.7/im and 7.9 /^m bands, leads to fluxes higher by up to 8% and 15% respectively. 
The smallest aperture is preferable in our case to minimize contamination by ghosts and other 
artifacts (see Section 3). It is also preferable, for the 3.6 /xm and 4.5 /im bands, because the 
underlying disk emission is non-negligible at these wavelengths (see Table 1). No extended 
emission is detected in the 5.7 and 7.9 /im bands, and total fluxes measured in a large aperture 
(with correction factors of 0.63 and 0.69, respectively, that are appropriate for extended 
emission) are lower than the nuclear fluxes by 25% and 6% respectively. Since the aperture 
correction factors for extended emission are much more uncertain than for point sources, we 
adopt the nuclear fluxes as total measurements in the 5.7 and 7.9 /im bands. For the total 
fluxes in the 3.6 and 4.5 /xm bands, we applied correction factors of 0.94 (IRAC handbook). 
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As the MIPS detectors do not suffer from the hght scattering phenomenon making accurate 
photometry difficult in IRAC images, ffuxes in the MIPS bands were simply measured by 
summation over the whole extent of the source. 



2.3. Optical imaging and spectroscopy 

Images in the B, V, R and I bandpasses were acquired at the CTIO 1.5 m telescope 
with the CFCCD (Cassegrain Focus CCD) instrument, in October 2001. The pixel size is 
0.43" and the angular resolution (full width at half maximum), measured on stars nearby 
NGC 1377, varies between 1.15" and 1.23" in the different bandpasses. For comparison 
between the B band and the I band, the images were registered using the positions of a 
dozen stars, with a residual radial dispersion of 0.2", and convolved to a common resolution 
of 1.5". 

Optical long-slit spectra between 3600 and 6900A, at a spectral resolution of SA, were 
acquired at the Steward Observatory 2.3 m telescope in November 2001. The observations 
and data reduction are described by Moustakas et al. (2006, in preparation). We used an 
extraction aperture of 2.5" x 10" centered on the nucleus of the galaxy and oriented along 
the minor axis. 



2.4. Flux calibration of the infrared spectra 

The default flux calibration of Spitzer spectral data is appropriate only for staring 
observations of point sources, in which the loss of flux due to the diffraction beam being larger 
than the slit width is indirectly corrected for, by matching staring observations of stars with 
spectrophotometric models. Additional steps are necessary to calibrate mapping observations 
of arbitrarily-shaped sources. Although NGC 1377 is a point source in the infrared (Sect. 3), 
we cannot use the default flux calibration because it was not perfectly centered in the slits. 
The flux loss fraction as a function of wavelength was derived from estimates of the beam 
proflle using the PSF simulator distributed by the Spitzer Science Center.^ The spectra 
of the galaxy were then corrected by this function to remove the added point-source flux 
falling outside the central slit, because this flux is already retrieved by mapping the source. 
For the low-resolution data, it is also necessary to correct for the fact that the default flux 
calibration is done by extracting stellar spectra within an aperture much smaller than the 



^http://ssc. spitzer. caltech.edu/archanaly/contributed/browse. html 



slit length. The conversion from e~/s units to flux density units using the full slit length 
was derived from highly redundant observations of the calibrating stars HR6688 (for the 
LL module) and HR7891 (for the SL module), mapped in the slit-parallel direction. The 
spectra of these stars, extracted in an aperture of the same width as that used for NGC 1377, 
were divided by model spectra from Decin (2004)^ to derive the flux calibration functions. 
These functions were adjusted with smoothly- varying polynomials before being used to divide 
the LL and SL spectra of NGC 1377. Broadband photometry in the IRAC, MIPS and 
ISOCAM filters was then simulated from the spectra, using the appropriate transmission 
curves, and compared with imaging photometry. We obtain excellent agreement at 24 /im 
and 5.7;um: F24(MIPS)/F24(LL) = 1.01 and F5.7(IRAC)/F5.7(SL) = 0.98. The agreement 
with ISOCAM data is also quite good despite uncertainties in the saturation correction 
in the images (for details, see Roussel et al. 2003): F6.7(ISOCAM)/F6.7(SL) = 0.88 and 
Fi5 (ISOCAM) /Fi5(SL - LL) = 0.97. Note that the spectral shape between 13 and 14.5 /xm 
is not well calibrated at low resolution, because of an artifact at the 10% level in the detector, 
which cannot be modelled at present. Comparison with SH data, which are not affected by 
this artifact, however shows good consistency. We also obtain F7 9(IRAC)/F7 9(SL) = 0.92. 
Since IRAC photometry is itself uncertain (see Sect. 2.2), this does not imply any significant 
error in the spectral flux calibration. 



3. Morphology and general properties 

NGC 1377 is a member of the unrelaxed Eridanus galaxy group (Willmer et al. 1989), at 
an estimated distance of 21 Mpc. In stellar light, it has the appearance of a regular lenticular 
galaxy, with a diameter of 1.8' (about llkpc) at the 25 mag/arcsec^ blue isophote, an axis 
ratio close to 2.0, and a position angle of 92° (de Vaucouleurs et al. 1991). Heisler & 
Vader (1994), based on optical imaging of NGC 1377 as part of a sample of galaxies with 
Feofmi > -^100 Atm; mentioned the presence of a dust lane along the southern part of the minor 
axis, in an otherwise featureless morphology. The far-infrared (40-120 /xm) luminosity of 
NGC 1377, as derived from the IRAS 60 and 100 /xm fluxes, is 4.3 x 10^ Lq. The center of 
the galaxy contains a compact molecular gas reservoir of 2 x 10^ M0 or more, about ten 
times more massive than the amount of gas expected to be associated with the dust seen in 
emission (see Section 6), if a Galactic gas to dust mass ratio is assumed (Roussel et al. 2003). 
The molecular gas properties are similar to those characterizing starbursts: the far-infrared 
to CO flux ratio is seven times higher than in normal galaxies and the gas is subthermal. 



^The data are made available at the web address: http://ssc.spitzer.caltech.edu/irs/calib/templ/ 
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Figure 1 shows the (B-I) color map derived from the optical images obtained for SINGS. 
The disturbance noted by Heisler & Vader (1994) is readily apparent, but seems more com- 
plex than a single dust lane. Two main orthogonal red features are surrounded by more 
diffuse red filaments. The color excess between the central feature and the normal light dis- 
tribution immediately to the north is approximately 0.9 mag. The residuals from the B-band 
image after subtraction of a smooth symmetric photometric model are shown in Figure 2. 
Some faint arcs of excess light are seen in the center, as well as to the north-east and to 
the south-east of the central regions. Negative residuals of large amplitude, that may be 
interpreted as regions of excess obscuration and thus excess column density of cold dust and 
gas, are seen in the nucleus and several pockets immediately to the south. In addition, a 
straight lane extends further to the south, and diffuse structures further to the west. These 
faint residuals are discernible out to about 2kpc in projected distance from the nucleus. 
The overall morphology suggests a recent merging event with one or several small bodies, in 
which the structure of NGC 1377 would have been preserved. As NGC 1377 is very deficient 
in atomic hydrogen gas, more robust evidence cannot be obtained from HI observations but 
would require high spectral resolution and high sensitivity observations of the stellar and 
molecular gas kinematics. 

The infrared images from 3.6 to 24 /im are shown in Figure 3. They are dominated 
by a very bright nucleus. At 3.6 and 4.5 /im, the underlying smooth emission from the 
old stellar disk is also visible, with the same distribution as the I-band light. The reticle- 
shaped artifact in the 4.5 /im image (with orthogonal positive and negative rows) is produced 
by the "multiplexer bleed" and "pull-down" effects, respectively, typically associated with 
bright point sources (IRAC handbook). At 5.7 /im, the stellar disk emission has fallen below 
the sensitivity limit, and only the nucleus is detected at longer wavelengths. It has the 
appearance of the point response functions of the IRAC and MIPS detectors, with bright 
diffraction spikes. The rows that are bright throughout the field of view of the 5.7 and 
7.9 /im images exemplify the "pull-up" and "optical banding" artifacts (IRAC handbook). 
Since existing corrections of these artifacts are only cosmetic, we did not apply them but 
masked out the affected rows and columns in our analysis. 

To decompose the nuclear 7.9 /xm source and determine whether its size is measurable, 
point spread functions (PSF) derived from observations of bright stars are not appropriate, 
because their cores are more extended than the observed core of NGC 1377, sampled at the 
native IRAC pixel size of 1.22". This is probably because the spectral slope of NGC 1377 
within the 7.9 /im filter bandpass is much steeper than that of stars (the flux density decreas- 
ing sharply from 7.8 /im to 9.5 /im, because of the deep flanking silicate absorption), and the 
PSF width depends on the spectral slope, increasing as the spectrum reddens. The model 
PSFs computed with the tool provided by the Spitzer Science Center, "S-Tiny Tim" , using 
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the spectrum of NGC 1377, are also inappropriate for our purpose, because they depend 
on pre-launch characteristics and are not accurate enough. Within the uncertainties on the 
PSF, the mid-infrared source in NGC 1377 is thus not resolved. Both this constraint and that 
obtained previously on the size of the 2.12/im molecular hydrogen line emission (Roussel et 
al. 2003) imply that the region of activity is less extended than ~ 1", i.e. ~ 100 pc. To look 
for asymmetries in the 7.9 /xm emission, we extracted one half of the map on the left side of 
the vertical symmetry axis of the diffraction and "pull-up" pattern (the tilted upper half in 
Figure 3), and subtracted the mirror image from the other half of the map. The result is 
shown in Figure 4. The residuals do not include any extended emission, but primary and 
secondary ghosts of the nucleus, as well as the real source labelled A. Given its detection in 
the optical bands and its approximate fluxes in the IRAC bands, this source contains both 
stellar and dust components, and is likely a distant galaxy. 



4. Optical spectroscopy 

Optical slit spectroscopy shows the [Nil] line at 6583 A and the [Sll] lines at 6716 and 
6731 A in emission without any detection of Ha (Kim et al. 1995). Our long-slit data confirm 
this result (Fig. 5). The [OlII], [Oil] and H/3 emission lines are not detected either. As the 
optical depth of the infrared source is extremely high (Sect. 6), we infer that [Nil] and [SlI] 
arise in the foreground. Even accounting for underlying stellar absorption in the Ha line, 
the [NII](6583A)/Hq; ratio is constrained to be larger than unity. Since the lines cannot be 
excited by an active nucleus, because they do not probe the nuclear regions, they indicate 
shock excitation (Baldwin et al. 1981). In addition, the low [SlI](673lA)/[SlI](6717A) ratio 
implies that they arise in regions of low electronic density (Dopita 1977). 

Nebular diagnostics of metallicity are not useable in NGC 1377. Stellar absorption 
indices of metallic lines and H/3, on the other hand, are measureable in the spectrum of 
NGC 1377. But they are degenerate with respect to age and metallicity, except for galaxies 
which have been passive for several gigayears, which appears not to be the case for NGC 1377: 
EW(Mgb) = 2.4A, EW(Fe527o) = 2.lA and EW(H/5) = 3.3A in the index definition of 
Worthey et al. (1994). However, the large quantity of CO gas and dust (Sect. 6), relative to 
the stellar emission, argues against NGC 1377 being a low-metallicity system. 

The optical spectrum, covering a large wavelength interval (between 3600 and 6900A), 
enables us to constrain the star formation history of NGC 1377. We used the population 
synthesis models of Bruzual & Chariot (2003) to fit the entire spectrum by a discrete sum 
of instantaneously- formed populations, as described by Moustakas & Kennicutt (2006), as- 
suming solar metallicity and a Salpeter (1955) initial mass function. The estimated stellar 
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extinction is Ay ~ 0.9 mag in the screen hypothesis. We find that the hght is dominated by 
populations older than 1 Gyr (~ 60 to 80% in the V band), and that the most recent star 
formation episode traceable in the optical occurred more than 500 Myr ago. This is fully 
consistent with the radio continuum upper limit indicating that the cosmic ray population 
generated by the last major star formation episode has decayed, implying that this episode 
occurred more than 100 Myr ago. 



5. Infrared spectral energy distribution and spectral features 

The entire infrared spectral energy distribution, including the low-resolution spectra, 
calibrated to reproduce the broadband photometry as explained in Section 2.4, is shown in 
Figure 6. To the far- infrared broadband fiuxes, color corrections between 1.00 and 1.08 were 
applied. It is immediately apparent that the infrared source is not only very hot and thus 
compact (in agreement with the fact that it is not spatially resolved), but is also strongly 
self-absorbed, as shown by the deep absorption bands from amorphous silicates at 10 and 
18 /im. The forbidden lines from HII regions and photodissociation regions that are bright in 
star-forming galaxies ([Sill], [SlII], [Nell], [Nelll]) are absent in NGC 1377 (Table 2), except 
for a very weak [Nell] line at 12.8 /xm, which will be discussed in Section 7. The energy peak 
lies between 30 and 50 /xm (Fig. 6b), which is extremely unusual, even for starburst galaxies. 
An even more extreme spectrum (with an energy peak around 15-20 /xm) has been observed 
in SBS 0335-052, a blue compact dwarf galaxy of very low metallicity (Houck et al. 2004b). 
Since SBS 0335-052 contains one of the most deeply embedded super star clusters known, it 
may be instructive to compare it with NGC 1377. Their spectra suggest that the geometry of 
the source may be similar in both systems, compact and sharply bounded. There are however 
important differences. The mid-infrared emission of SBS 0335-052 between 5 and 10 /im is 
very faint, indicating that dust species related to the aromatic band carriers are absent and 
that very small grains are depleted with respect to big grains (Plante & Sauvage 2002). On 
the contrary, the mid- infrared emission of NGC 1377 accounts for a large fraction of the 
total infrared power and implies that very small grains are abundant (even though aromatic 
band carriers are seemingly absent), or alternatively that the radiation field intensity is high 
enough to heat big grains to very high temperatures. The optical depth in the amorphous 
silicate absorption bands is also much higher in NGC 1377. 

No other spectral features are seen, except some molecular hydrogen rotational lines 
(Fig. 7), which will be discussed in Section 8, as well as a few broad emission features of 
very low intensity and weak absorption features from gas-phase molecules. Superposed on 
the deep 10 /xm silicate absorption band, a very weak 11.3 /xm emission band is detected, and 
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likely arises from polycyclic aromatic hydrocarbons (Fig. 8). This feature is well matched by 
the profile of the 11.3 /im PAH band of the Orion bar (Peeters et al. 2004), at the spectral 
resolution of our data. It may originate from intermediate dust shells, with optical depths 
lower than the silicate optical depth, but still higher than unity (see below). The 11.3 /xm 
aromatic band is always observed in association with other bands, the most prominent of 
which lies at 7.7 /xm. The non-detection of the 7.7 /xm band could easily be due to its being 
hidden in the observed broad 5-9 /xm bump, which peaks almost exactly at the central 
wavelength of the aromatic band (see Section 6). Likewise, a 12.7/im aromatic band, if 
present, could be hidden in the broad bump on the red shoulder of the 10 /xm silicate band 
(Fig. 6). 

Another broad and faint feature is seen below the 17.03 ;um H2 line, at the bottom of the 
18 /im silicate absorption band (Figure 8). It could be composed of the same bands at 17.1, 
17.4 and 17.8 /im that were shown recently to be common in the interstellar medium of normal 
galaxies (Werner et al. 2004; Smith et al. 2004, 2006), and that are thought to be related to 
polycyclic aromatic hydrocarbons. Van Kerckhoven et al. (2000) have extracted a plateau of 
blended broad features attributed to the C-C-C bending modes of PAHs in this wavelength 
range (more exactly between 15 and 20 /im), in compact HII regions as well as young and 
evolved stars. They also state that small-size PAHs would have sharper spectral structure 
(in particular a band at 16.4 /xm, unseen in NGC 1377) than bigger aromatic compounds. If 
the 16.5-18.5 /im bump observed at the bottom of the 18 /im silicate band in NGC 1377 is 
attributable to PAHs, then the lack of substructure would indicate a predominance of bigger 
aromatic compounds. 

Finally, we tentatively detect two additional broad features at 6.3 and 6.7/im (Fig. 8). 
The dip longward of 6.8 /im could be caused by an absorption from C-H bending modes of 
hydrocarbons at 6.85 /im, seen in deeply obscured infrared-luminous galaxies by Spoon et 
al. (2001) and Spoon et al. (2004). For comparison. Fig. 8 shows the spectrum of one of 
these galaxies, IRAS 08572-1-3915 (Spoon et al. 2006a), where the absorption bands from 
hydrocarbons at 6.85 /im and 7.25 /im are clearly seen. However, there is no known species 
being able to cause the dip at 6.5 /im. The 6.3 /im band has been observed in at least two 
post-AGB stars by Peeters et al. (2002) and has also been invoked by Sturm et al. (2000) to 
explain the excess emission redward of the 6.2 /xm aromatic band of starburst galaxies. It has 
so far not been observed alone (without the 6.2 /^m band) in galaxies. According to Peeters 
et al. (2002), the class of PAHs with a 6.3 /im peak wavelength may have a more pristine 
composition than the 6.2 /im PAHs, and may not yet have been exposed to hard radiation. 
It can be expected from the compact starburst nature of NGC 1377 that PAHs have all been 
destroyed in the zone of activity and survive only in outer shells where they would not yet 
have felt the influence of the starburst ; it would however be difficult to understand how 
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they could have avoided some kind of processing. The as yet unidentified 6.7;um band has 
been observed in NGC7023 (Werner et al. 2004) and in spiral galaxies (Smith et al. 2006). 

If we assume that the intrinsic relative intensities of the 6.3, 6.7, 17-18 and 11.3 /xm 
bands are similar to what is seen in normal galaxies, then they are not consistent with zero 
extinction, but would be consistent with an optical depth of the order of half the optical 
depth affecting the continuum (which is about 8.5 at 11.3 /xm, as constrained in Sect. 6). 
Given the utter uncertainty on the relative geometry of the different dust phases, and the 
fact that the relative intensities of the aromatic bands are intrinsically variable in various 
interstellar sources and in galaxies (Peeters et al. 2002; Smith et al. 2006), the identification 
of the reported features with PAHs is therefore plausible. 

Very faint absorption features are seen at high resolution between the two amorphous 
silicate bands: HON at 14.03 /im, which lies at the end of two spectral orders but is confirmed 
in both, possibly C2H2 at 13.7/im, and CO2 at 15.0 /im, which lies at the end of a single 
spectral order and is outside the wavelength range of the adjacent order (Fig. 9). These 
absorption bands, tracers of hot dense gas, have been observed in deeply obscured infrared- 
luminous galaxies (Spoon et al. 2006b). Contrary to what is seen in some of the latter objects 
and in Galactic center sources (Chiar et al. 2000), we do not detect any absorption features 
in the 5-8 /im range, i.e. from water ice at 6/im and from hydrocarbons at 6.85 and 7.25 ;um. 



6. Geometry and density 

To attempt to constrain the geometry of the infrared source in NGC 1377, following the 
same approach as Plante & Sauvage (2002) for SBS 0335-052, we used the DUSTY radiative 
transfer model (Nenkova et al. 2000), assuming a unique, spherically-symmetric source. By 
construction, aromatic compounds are not included in the model, and impulsive heating^ 
of very small grains is not taken into account. This shortcoming is partially mitigated by 
the fact that the dust in NGC 1377 is heated to very high temperatures, and thus closer to 
thermal equilibrium than in normal starbursts, but the results should be interpreted with 
caution. However, we note that the effect of neglecting impulsive heating should produce 
residuals smoothly distributed in wavelength, not narrow features (see below). Another 
important limitation is that the same grain size distribution has to be postulated for all 
dust species, despite the fact that all the families of dust models of the local interstellar 
medium need very different size distributions for silicates and carbon-based dust (Zubko et 
al. 2004). Since no unique set of geometric parameters will provide a fit to the data, and 



^also referred to as transient, stochastic, or single-photon heating 
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extensive exploration of the parameter space would be prohibitively long for such a high 
optical depth as seen in NGC 1377, we restrict the analysis to a qualitative assessment of the 
density. For the silicate grains, we find that, among the choices available within DUSTY, the 
optical properties of Ossenkopff et al. (1992) allow the best reproduction of the shape of the 
silicate absorption bands. For the carbonaceous grains, graphite, with the optical properties 
of Draine & Lee (1984), allows a better fit to the spectrum of NGC 1377 than amorphous 
carbon. 

Figure 10 shows fits reproducing the available data from 3.6 to 156 /xm, except in a 
few wavelength intervals that will be discussed separately below. The models contain at 
least 70% of amorphous silicates by number of grains, which corresponds to about 80% of 
silicates by mass, using the mass densities provided by Weingartner & Draine (2001) and 
Laor & Draine (1993). This mass abundance is a bit larger than required by the relevant 
class of dust models of Zubko et al. (2004) (~ 66%), as well as by the dust model of Li & 
Draine (2001) (~ 73%). Such a large fraction of silicates is made necessary in NGC 1377 
to account for the great depth of the 10 /xm and 18 /im features. It can be speculated that 
this refiects composition changes effected by the intense radiation in the starburst region, 
silicates possibly being more resilient than carbonaceous grains in an oxygen-rich medium. 
The best fits are also obtained for a grain size distribution enhanced in small grains, with 
a power-law index —5 instead of —3.5 (for the standard Mathis et al. (1977) distribution) 
and a lower cutoff of the order of 10 A. In particular, for a normal size distribution, the 
opacity would have to be much higher in order to reproduce correctly the 10 /im and 18;Um 
absorption bands. This may be an artifact of the model, compensating for the absence of 
impulsive heating. 

The model implies an opacity of at least 75 in the V band, more than 20 at the center of 
the 10 /im silicate band, and about 10 in the 18 /im silicate band. The radial density profile, 
that was parameterized as a broken power-law, is very steep, with an index of —2 in the 
center and —1.5 outside. The density profile in the outer parts is ill-constrained, and does 
not allow us to derive a total size. For a bolometric luminosity of 1.1 x 10^'^ Lq, derived from 
the data themselves by integration over the whole infrared spectral energy distribution, and 
a temperature of 2000 K at the inner surface of the dust shell, the radius of this inner layer is 
constrained to be 0.08 pc, very similar to what Plante & Sauvage (2002) found for SBS0335- 
052, although the total infrared luminosity of NGC 1377 is about three times higher (the 
inner dust temperature of the model of SBS 0335-052 is only 700 K). The dust mass derived 
from the 60 /im and 100 /im fiuxes, using the Galactic emissivity as formulated by Bianchi 
et al. (1999), is 2.7 x 10^ Mq. It can alternatively be derived from the opacity obtained with 
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the DUSTY model: 

TV = ^(/i fi^Vi Pi) / (4/37ra^ n{a) da) / {N{r) dr) 

■ J a J r 

and Mdust = ^U\ Pi) / (4/3™^ n{a) da) {N{r) ATrr^dr), 

■ J a J r 

where the index i designates each dust species (here sihcates and graphite), f\ is their fraction 
by number, /tyi is their mass absorption coefficient in the V band, pi is the sohd density 
of the grains, n{a) da the grain size distribution and N{r) dr the radial density distribution 
adopted in the model. It follows that 

M 4^. Ei(/iPi) L{N{r)r^dr) 

Eil/i/^viPi) J^{N{r)dr) 

We used /ty = 0.3m^g~^ for silicates and 5m^g~^ for graphite (Draine & Lee 1984) ; and 
p = 3.5gcm~'^ for silicates and 2.24 gcm~^ for graphite (Weingartner & Draine 2001; Laor 
& Draine 1993). Varying the V-band optical depth between 75 and 95 (Fig. 10), we obtain a 
dust mass in the range (2.6-3.3) x 10^ M©, in good agreement with the determination from 
the far-infrared fluxes. Assuming a gas to dust mass ratio of 100 (Sodroski et al. 1994), we 
obtain a central gas density of about 3 x IO^MqPc"^, or equivalently a hydrogen nucleus 
density of lO^cm"^. 

It should be noted here that the observations up to 156 /im are not sensitive to the 
potential presence of very cold dust (not illuminated by the starburst). Our dust mass 
estimates are therefore lower limits, and we suspect the total dust mass to be about ten 
times higher, since the mass of cold molecular hydrogen is ten times higher than the mass 
of gas expected to be associated with the dust seen in emission. 

The residuals, after subtraction of the model (Fig. 11), can tentatively be identified with 
emission bands from aromatic compounds and absorption bands from crystalline silicates. 
If the identification with aromatic band carriers is correct, their chemical composition is 
likely different from those dominating the mid-infrared spectrum of normal star-forming 
galaxies, as indicated by the absence of the 6.2 /im band near the 6.3 /im band (Section 5). 
In Figure 11a, prominent features are seen at peak rest wavelengths of 7.7 and 12.7pm. Note 
that the absorption profile of silicates is not perfectly reproduced ; it might be possible to 
explain the excess emission by a combination of impulsive heating (neglected in the model) 
and modified optical properties of silicates. 

Absorption bands from crystalline silicates have been found in some infrared-luminous 
galaxies by Spoon et al. (2006a). In Figure lib, features are seen at rest wavelengths of 
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23.5 and 28/im. Other known features from forsterite exist at 16 and 19.5 /im, but the 
quantification of their depth in NGC 1377 is made hazardous by the possible contamination 
by aromatic bands in the 17-18 /im range (Fig. 8). The 16 /im band, in particular, would 
have to be shifted to longer wavelengths in comparison with IRAS 08572+3915. Another 
feature exists at ~ 11 /im, but there again it is confused by the 11.3 /im PAH. Crystalline 
forsterite also produces strong bands at about 34 and 70 /^m, among others (Koike et al. 
1993; Molster et al. 2002). The dip of the spectrum of NGC 1377 at the long-wavelength end 
of the LLl module, between 30 and 35 /im (Fig. 6 and 11) may be attributed to absorption 
by crystalline silicates. Likewise, they may contribute, at least partially, to the weakness of 
the 71 //m broadband fiux relative to the 60, 100 and 156 //m fiuxes. 



7. Forbidden lines 

Mid-infrared forbidden lines are important tracers of HII and photodissociation regions, 
and of the physical conditions therein. Figure 9 shows the detection of a very weak [Nell] line. 
[Nell] emission without detection of any other forbidden line in the mid-infrared range is also 
observed in some galaxies of the sample of Spoon et al. (2006a). We attempt here to extract, 
from the detection of [Nell] and the absence of other lines, constraints on the recombination 
rate and electronic density. The low signal to noise ratio of [Nell] precludes measuring the 
line width accurately, but it is well matched by the spectral resolution (corresponding to a 
full width at half maximum of 500kms~^ for an unresolved line). 

The estimated fiux (Table 2) can be translated into a recombination rate. We have 

F[Neii] X AixD^ = \iuV C Ar(Ne+) N{e-) {Tjl K)-^-^ n{T,) g^' exp{-hu / (kT,)) 

and Nr = a(Tc) V N(R+) iV(e~), where D is the distance of NGC 1377, V is the volume of the 
HII region, C = 8.6287 x 10~^cm^s~^, A^ designates a volume density, f2(Tc) is the collision 
strength, g\ is the statistical weight of the ground level (4 for [Nell] and 5 for [Nelll]), and 
a(Te), which depends only weakly on the electronic density, is the hydrogen recombination 
coefficient. We assume that Ar(Ne)/A^(H) = [Ar(Ne+) + Ar(Ne2+)]/Ar(H+) and that the 
neon abundance is equal to the solar value, 6.9 x 10^^ (Asplund et al. 2004). The collision 
strengths of the [Nell] and [Nelll] transitions are tabulated respectively by Saraph & Tully 
(1994) and Butler & Zeippen (1994), and the hydrogen recombination coefficient by Storey 
& Hummer (1995). Accounting for the upper limit on the [Nelll] line fiux and assuming that 
the electronic temperature is between 5 x 10^ K and 10^ K, we obtain a recombination rate 
(uncorrected for the extinction in the [Nell] and [Nelll] lines) N,- = (5.7 ± 1.6) x 10^^ s~^. 

The upper limit on the thermal emission at 3 cm (Roussel et al. 2003) independently 
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provides an upper limit on the recombination rate: N,- < (Tc/10^K)^°'^^ x 5.4 x 10^^ s~^, 
which is at least an order of magnitude smaller than the intrinsic ionizing photon flux A^Lyc, 
because most of the ionizing photons are absorbed by dust. The estimate derived from 
the [Nell] line is consistent with this limit provided that the optical depth in the line be 
lower than 2.6 . The modelled optical depth at the wavelengths of both [Nell] and [Nelll] 
is however of the order of 6.5 (Sect. 6). This discrepancy calls for at least one remark: if 
the starburst is not made of a unique stellar cluster but occurs simultaneously at several 
locations (Sect. 6), then the [Nell] line emission may arise from the regions of lowest optical 
depth. 

The upper limits obtained on the hydrogen recombination lines (Roussel et al. 2003) 
provide another constraint on the nebular extinction. The requirement that the recombina- 
tion rates derived from the 3 cm, [Nell] , Br7 and Pa/? measurements be all consistent implies 
a nebular optical depth of 0.9 < r(12.8/im) < 2.6, assuming the Galactic center extinction 
law of Moneti et al. (2001). 

We have assumed that collisional de-excitation is negligible. If this were not the case, 
then our estimate of the apparent ionizing photon flux would be underestimated, and the 
optical depth in the [Nell] line would have to be even lower. The data do not provide any 
constraint on the electronic density, but the radiative transfer model (Sect. 6) suggests a 
central gas density of the order of lO^cm"^. To determine whether the mid-infrared nebular 
lines can be suppressed by collisional de-excitation in NGC 1377, we computed the critical 
densities of the ionic lines that are usually the brightest between 10 and 35 /xm in star-forming 
galaxies, and that are tabulated in Table 2. If the electronic density is below the critical 
densities, then we expect minimum line ratios -Fpiii] /-^[NcII] ~ 0.5, -F[siii]33.5^m/-^[Neii] ~ 0.15 
and -F[siii]i8.7^m/-^[Ncii] ~ 0.2. These values are derived empirically from the lower bounds 
observed in nuclei within the rest of the SINGS galaxy sample ; they correspond to the 
theoretical values if ~ 60% of gas-phase silicon is in Si^ assuming that ^ 90% of silicon 
atoms are retained in dust grains, and ^ 5% of sulfur is in S^^, with the Si and S solar 
abundances given by Grevesse & Sauval (1998). The upper limits obtained for NGC 1377 
are above these minimum ratios, except -Fpm] i8.7^m/-^[Neii] < 0.12. Since an optical depth of 
~ 1 in [Nell] (corresponding to an optical depth of the same order in [Nelll]) is sufficient to 
explain this low ratio, we are unable to set constraints on the electronic density. Collisional 
de-excitation may not operate for [Nell] and [Nelll], since their critical densities are at least 
ten times higher than the neutral density estimate from the radiative transfer model. It 
should be noted, on one hand, that the peak density could be significantly higher than our 
estimate, and, on the other hand, that the electronic density should be much lower than 
the neutral density if pressure equilibrium is maintained. Is is thus unlikely that [Nell] and 
[Nelll] be collisionally de-excited. 
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8. Molecular hydrogen 

NGC 1377 is the only member of the SINGS sample in which H2 lines dominate the line 
spectrum in the mid-infrared (and also in the near-infrared) . Three rotational transitions of 
H2 have been confidently detected: the S(l) and S(2) transitions at high spectral resolution, 
and the S(3) transition at low resolution in the SLl module (Fig. 7). The line fluxes are 
given in Table 2. Based on previous observations of rovibrational transitions in the near- 
infrared, we had concluded that H2 is collisionally excited (pure fluorescence being ruled out 
by the high flux ratio of the v=l-0 S(l) transition to the v=2-l S(l) transition), and that 
slow shocks in the molecular phase are responsible (Roussel et al. 2003). The alternative 
mechanism, collisions with hot hydrogen atoms and molecules in dense photodissociation 
regions (Hollenbach & Tielens 1997, and references therein), seems indeed in conflict with 
the expectation that photoelectric heating will be inefficient in NGC 1377, because small 
dust grains may be positively charged due to the very intense radiation field. This effect 
was observed by Malhotra et al. (2001) for the [Cll] and [Ol] lines at 157.7;um and 63.2 /im, 
which are among the most efficient coolants of photodissociation regions. They were not 
detected by ISO-LWS observations in NGC 1377 (Roussel et al. 2003). The power emitted 
in the sum of the S(l) to S(3) transitions of H2 represents a fraction ~ 6 x 10^'^ of the 
far-infrared luminosity (40-120 /im) ; and the upper limits on the [Cll] and [Ol] luminosities 
respectively ~ 6 x 10~^ and ~ 9 x 10~^. 

Using the recent photodissociation region models of Kaufman et al. (2006), we can test 
whether this excitation mechanism is efficient enough to account for the observed bright- 
ness of the H2 lines. In these models, the maximum fraction of the far-ultraviolet radi- 
ation (FUV, between 6 and 13.6 eV) emerging in the rotational S(l) line is ~ 10~^, and 
it occurs for a radiation field intensity of Gq ~ 200 times the local value and densities of 
n ~ lO'^cm"^. At higher Go values, the gas heating efficiency declines ; and at higher densi- 
ties, the gas cools partially through collisions with dust grains (Hollenbach & Tielens 1997). 
The H2 line flux ratios of NGC 1377 (Fs(i)/Fs(o) > 7, Fs(2)/Fs(o) > 4, Fs(3)/Fs(i) > 1 and 
Fv=i_os(i)/-^v=2-is(i) > 30, these limits accounting for any amount of extinction) indicate 
that Go > 10"^ and 10^ < (n/cm^^) < 10^ ; and under these conditions, the maximum 
Fs{i)/ FUV flux ratio is ~ 4 x 10~^. In NGC 1377, the 3cm upper limit implies that the 
part of the ionizing photon flux A^Lyc escaping dust absorption is < 7.4 x 10^^ s~^, and 
Lp\jy ~ 5.1 X 10~"'^^W X {Ni^yc/s~^) for a 1 Myr-old starburst with a Salpeter initial mass 
function (using the population synthesis model of Leitherer et al. (1999)). Therefore, we 
expect to observe a maximum of ~ 30 x 10^^^ Wm^^ in the v=0 S(l) line. 

The observed flux is twice as high. The intrinsic flux can be much higher, if it is affected 
by the same optical depth as toward [Nell], since gas heated in photodissociation regions 
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would have to be located very close to the HII regions. The minimum optical depth toward 
[Nell] (Section 7) would correspond to r(17/im) ~ 1.1. The intrinsic flux in the H2 v=0 S(l) 
line would therefore be about six times higher than the model prediction. This difference is 
marginal, compared with model uncertainties, in particular on the H2 coUisional de-excitation 
rates. However, we note that the Br7 to H2 v=l-0 S(l) flux ratio (which is immune to 
extinction effects if the emission is nearly cospatial) is lower than 0.1 in NGC 1377 (Roussel 
et al. 2003), which is much lower than expected in the case of photodissociation region 
heating of H2. For example, the normal starburst NGC 1022 has a Br7 to H2 v=l-0 S(l) 
flux ratio of 2.2 . We thus argue that shock excitation plays a dominant role in NGC 1377, 
although photodissociation region models are marginally consistent with the present data. 

An excitation diagram combining all the observed transitions to date (Fig. 12) conflrms 
collisional excitation, and allows us to set a constraint on the amount of extinction in the 
H2 lines, as explained below. The flux of a transition f^f/J^J — 2 is given by 
F = hv A A^ucoi ^/(47r), where N^^oi is the column density of molecules in the upper level, A 
is the spontaneous emission probability, hv is the transition energy and Vt is the source solid 
angle. We also have A^'^ = g^ Ntot exp(-Eu / (kT)) / Z{T), where g^ = (2i + 1) (2J + 1) is 
the statistical weight (with the spin number i = for even J or para transitions, and i = 1 
for odd J or ortho transitions), Z(T) ~ 0.0247 T / (1 — ea:p(— 6000K / T)) is the partition 
function, and E^ is the upper level energy (Herbst et al. 1996). Line ratios thus allow the 
derivation of an apparent excitation temperature: 

kT = (E„2 - ^ui) / ln(iVui/A^u2 X gu2/9ui) 
with N^i/N^2 = Fi/F2 X A2M1 X A1/A2. 

We used the spontaneous emission probabilities given by Turner et al. (1977). Since the 
temperature constrained by the rotation- vibration S(0) lines and by the rotation- vibration 
S(l) lines is consistent with the temperature derived from the pure-rotation transitions, 
collisional excitation is conflrmed. 

In order for the excitation temperature to remain monotonic as a function of the upper 
level energy (T increasing as E^^ increases), the equivalent V-band optical depth has to be 
lower than 17 (r < 1.1 at 17/im). This estimate takes into account the errors on the line 
fluxes, and assumes that H2 at all temperatures is affected by a same extinction following 
the Galactic center extinction law of Moneti et al. (2001). Since the optical depth in the dust 
phase is much higher, most of the detected molecular hydrogen must come from intermediate 
layers. For zero extinction, the excitation temperature varies from (365 ± 50) K for the low- 
energy transitions, to (610 ± 25) K between the v=0 S(3) and v=l-0 S(l) transitions, to 
< 1500 K for the higher-energy transitions ; for the maximum allowed extinction, it varies 
between 420 K and < 1350 K. 
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Coinputing the mass of warm molecular hydrogen from the flux of the S(l) transition and 
the temperature derived from the rotational lines, we obtain 7.2 x 10^ Mq in the hypothesis 
of zero extinction, and 1.9 x 10^ Mq with ry = 17. Since we had derived a total molecular 
hydrogen mass of the order of 2 x 10^ Mq from CO observations (Roussel et al. 2003), this 
implies that less than 1% of H2 is heated to more than 300 K. In comparison with the 
starbursts observed by Rigopoulou et al. (2002), this fraction is low. However, in the latter 
objects, the contribution of photodissociation regions to the heating of H2 is likely dominant, 
and the volume fraction of the molecular gas directly exposed to FUV radiation is probably 
much larger than in NGC 1377. 



9. Alternative hypotheses to a nascent starburst 

Different scenarios to account for the synchrotron deficiency (starburst deficient in high- 
mass stars, abnormally weak or strong magnetic field, dominant inverse Compton losses, 
extremely high electron opacity) were discussed at length by Roussel et al. (2003), and it was 
concluded that the most likely explanation involved a very young starburst. The possibility 
that a substantial population of relativistic electrons exists, but decays at a higher rate 
than the type-II supernova rate, in the hypothesis of inverse Compton losses dominating 
synchrotron losses, is not formally ruled out, but seems in conflict with other properties 
of the starburst. In particular, we did not detect the [Fell] line at 1.644 /xm, which is a 
quantitative tracer of shocks in supernova remnants (Vanzi & Rieke 1997), where cosmic 
ray electrons are accelerated ; the line flux ratio of [Fell] to II2 at 2.122/im is below 0.1 
in NGC 1377, whereas it is above 2 for the mature starburst in NGC 1022, which we chose 
as a comparison galaxy (see also Larkin et al. (1998) and Dale et al. (2004)). Reconciling 
these ratios by invoking differential extinction effects (and assuming cospatial emission) 
would imply ry > 50, whereas we have derived ry < 17 for H2 (Section 8). The massive 
stars in NGC 1377 are thus unlikely to have reached the supernova stage. In addition, it is 
hardly conceivable that a starburst older than 5 Myr would still be completely embedded, 
despite the effects of winds from massive stars and supernova shock waves on the interstellar 
medium. On the other hand, most cosmic rays generated by past star formation episodes 
would presumably have propagated over large distances and be located outside the current 
starburst region, which is very compact, and thus outside of reach of the intense radiation 
field of the starburst. Their absence is thus more plausibly due to NGC 1377 being previouly 
quiescent for a long period, than to their suppression by inverse Compton losses. Finally, 
we know of some ultraluminous galaxies, studied by Condon et al. (1991), that host very 
compact starbursts with radiation fields at least as intense as in NGC 1377 ; yet, they are 
bright synchrotron emitters, which implies that their magnetic field has been amplified to 
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the degree necessary to balance inverse Compton losses by synchrotron losses. It is unclear 
why magnetic field amplification would be inhibited in NGC 1377. In summary, the issue of 
inverse Compton losses is not settled by the present data, but we favor the hypothesis of a 
nascent starburst as much more likely. 

Simple energetics arguments also led Roussel et al. (2003) to conclude that it would 
be exceedingly difficult to account for the observed properties of NGC 1377 by accretion 
onto a massive black hole. An upper limit on the mass of a black hole at the center of 
NGC 1377, of the order of 2 x 10^ Mq, is derived consistently both from the lack of radio 
continuum emission (using an empirical relation between the minimum radio power emitted 
and the black hole mass of radio-quiet active nuclei) and from the narrow velocity profile 
of the H2 2.12 /im line (assuming virial equilibrium). In the unrealistic assumption that the 
efficiency of dust heating by the accretion power be 100%, this power would still have to 
be above the Eddington limit for a black hole of such a mass, to account for the observed 
infrared luminosity. In addition, substantial free-free emission would be observed, because a 
large fraction of the bolometric luminosity of active nuclei is ionizing ; the expected fiux, at 
wavelengths short enough that thermal opacity is negligible, is well above the observational 
upper limit. A new element is brought by the Spitzer data to constrain the nature of the 
power source. The fact that no high-excitation lines are detected makes the hypothesis of 
non-stellar activity all the more difficult. In particular, since the [Nell] line is detected, 
the undetected [Nelll] line, affected by almost the same extinction and having a similar 
critical density, sets a meaningful constraint on the radiation hardness. Based on the results 
of Sturm et al. (2002), pure Seyfert nuclei (i.e. minimally contaminated by star formation) 
have -F[Nciii]/-^[Ncii] ratios between 1 and 3. This line flux ratio is lower than 0.5 in NGC 1377. 
We conclude that non-stellar activity is very unlikely to be responsible for the infrared power 
of NGC 1377. 



10. Summary and discussion 

NGC 1377 first attracted our attention as being extremely deficient in radio continuum 
emission, with respect to its infrared brightness. Its dust is also heated to very high temper- 
atures, which implies an intense radiation field. The lack of synchrotron emission points to 
the galaxy being previously quiescent for at least 100 Myr (which is corroborated by optical 
spectral synthesis in Section 4), and the current star-forming episode being younger than the 
lifetime of type-II supernova progenitors (Roussel et al. 2003). Its far-infrared (40-120 /xm) 
luminosity is about 4 x 10^ L©, and the bolometric luminosity of the starburst is of the 
order of 1.1 x 10^° Lq. This corresponds to a starburst with a 1.6 x 10^ M© stellar mass. 
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assuming an age of 1 Myr and a Salpeter initial mass function between 0.1 and 120 M0 (or 
6 X 10^ Mq for a lower mass cutoff of IMq), and about 3 x 10'^ O stars. The global cold 
molecular gas properties are similar to those of starbursts, with subthermal excitation and 
a high far-infrared to CO flux ratio. The molecular gas mass, ~ 2 x 10* M©, is only 10 to 
30 times higher than the estimated starburst mass, which implies an extremely short gas 
exhaustion timescale for NGC 1377. The size of the infrared source is of the order of 100 pc 
or smaller, and the new observations reveal that the optical depth in the amorphous silicate 
absorption bands is one of the highest observed in extragalactic systems: the apparent opti- 
cal depth at 10 /im, of the order of 4, is very similar to that of IRAS 08572-1-3915, the most 
deeply obscured object in the sample of infrared-luminous galaxies of Spoon et al. (2006a). 
It should be noted that the modelled optical depth of NGC 1377 is about five times higher 
than the apparent depth. A few absorption bands from dense and warm molecular gas are 
also present, although very faint. 

The new data confirm that NGC 1377 has the characteristics of an extremely young, 
totally opaque nuclear starburst. It is in such an early phase of this starburst that the 
only tracer of HII regions yet detected is a very weak [Nell] line. Due to a combination 
of two effects, suppression of the ionizing photons by dust and high optical depth in the 
mid-infrared, the [Nell] line probes only ~ 0.8% of the intrinsic ionizing photon flux of the 
starburst. The growth of the HII regions is probably still inhibited by the pressure of the 
dense interstellar medium, and a large fraction of the massive stars must be still embedded 
in their parent molecular cloud. In order to reconcile the presence of very massive stars in 
nuclear starbursts with their nebular line ratios typically indicating low excitation (Thornley 
et al. 2000), Rigby & Rieke (2004) proposed that in such environments, high pressure causes 
massive stars to spend a large fraction of their lifetime in ultracompact HII regions, where 
the production of nebular lines is quelled. NGC 1377 could represent an extreme example 
of this phenomenon. Galactic systems where the dust competes so efficiently with gas for 
absorption of the stellar radiation are extremely rare. 

If we try to understand the starburst in NGC 1377 as a scaled-up version of stellar- 
size objects (as a mental picture only), the class of Becklin-Neugebauer (BN) objects may 
represent the closest analogs suited to that purpose. BN objects are at the transition between 
massive protostars and pre-main-sequence stars exciting ultracompact HII regions, after the 
intense accretion phase has stopped, but before the radiation pressure has balanced the 
interstellar medium pressure (Henning & Giirtler 1986). The emergent infrared radiation is 
strongly self-absorbed like in NGC 1377, and molecular hydrogen is excited by shocks, caused 
by the interaction of stellar outflows with the molecular envelope. It is unclear whether the 
same H2 heating mechanism applies to NGC 1377. Shocks could be caused by the collective 
outflows and winds of massive protostars and young stars ; or they could be caused by a 
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merger event, provided the disturbances revealed by optical color maps (Sect. 3) may be 
interpreted in this fashion. The optical forbidden line spectrum supports the existence of 
shocks also in the foreground of the region of activity (Sect. 4). 

The central neutral density inferred from a radiative transfer model is of the order of 10^ 
hydrogen nuclei cm~^, which is not greater than typical densities of giant molecular clouds. 
The estimated starburst mass (see above) is also very similar to the mass of the super star 
cluster coincident with the brightest mid-infrared source in the Antennae galaxies (Gilbert 
et al. 2000), although the latter cluster is more evolved, since it is associated with bright 
HII and photodissociation regions. In view of these results, we are not able to determine 
whether the activity in NGC 1377 arises from the formation of a single cluster, or from several 
isolated pockets. We had been led to the latter hypothesis by the fact that the molecular 
gas could be unstable with respect to gravitational instabilities (Roussel et al. 2003), so 
that star formation could be triggered almost instantaneously at several locations. The data 
are however compatible with the formation of a single stellar cluster in the nucleus. This 
question should be more easily addressed for more luminous objects. 

At the other end of the mass and luminosity scale, the infrared spectrum of NGC 1377 
is indeed strikingly similar to those of some deeply obscured infrared- luminous galaxies 
(Fig. 13). NGC4418 and IRAS 08572+3915 for instance (Spoon et al. 2001, 2006a) share 
the main spectral characteristics of NGC 1377. In addition, these luminous galaxies are 
part of a statistical sample of nascent starbursts that we have selected by their deviation 
from the infrared-radio correlation and their hot dust. Forthcoming papers will present and 
discuss the infrared-radio properties and molecular gas content of this sample, based on 
observations with Spitzer, the VLA array and the IRAM 30 m telescope. The characteristics 
of these galaxies will be contrasted with those of NGC 1377, serving as an archetype of the 
nascent starburst class. Their diverse radio continuum properties suggest that they should 
sample well young starbursts at various stages of development (Roussel et al. 2005). It will 
be one of our aims to sketch an empirical evolutionary sequence from dust-bounded nascent 
starbursts such as NGC 1377 to mature starbursts such as M82. 

We have previously derived a rate of occurrence of the nascent starburst phenomenon of 
~ 16% among galaxies whose flux density distribution peaks at 60 /im (Roussel et al. 2003). 
This rate is not in conflict with the expected timescales of star formation. Nascent starbursts, 
as we have deflned them, must fulflU two conditions: their current star formation episode 
must be younger than 5 Myr, and they must previously have been quiescent for ~ 100 Myr. 
If we note the typical duration of compact starbursts t x 5 Myr, where t is probably not much 
larger than unity, then the fraction of starbursts occurring in isolated mode (as opposed to 
recurrent mode) would be 0.16 x t. 
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Starburst tracks currently do not incorporate the pre-HlI nucleus phase. Likewise, the 
sets of galactic templates used to interpret the infrared counts in deep surveys do not include 
objects such as NGC 1377. Since even small changes in the assumed spectrum of aromatic 
bands in the mid-infrared yield important differences in the interpretation of the counts, as 
emphasized by Smith et al. (2006), the consequences of the existence of luminous galaxies 
sharing the properties of NGC 1377 cannot be neglected. This is illustrated by Figure 14, 
showing the large differences in the redshift-luminosity curve for the 24 /im MIPS filter, ac- 
cording to the assumed mid-infrared spectral energy distribution: the curve obtained for 
NGC 1377 is vastly different from those of popular starburst templates. The nascent star- 
bursts (selected by their deficiency in radio continuum emission, with respect to the infrared) 
are rare objects, in agreement with the brevity of their evolutionary phase. However, they 
represent a non-negligible fraction of galaxies with Fgo^m > -^loo^mj whose frequency may 
be much higher at earlier epochs, as suggested by the results of Chapman et al. (2004). 

We wish to thank Marc Sauvage for fruitful discussions, Mark Wolfire for tirelessly 
running photodissociation region models to help interpret the data, and Nancy Silbermann 
for the timely scheduling of the Spitzer IRS observations. Support for this work, part of the 
Spitzer Space Telescope Legacy Science Program, was provided by NASA through an award 
issued by the Jet Propulsion Laboratory, California Institute of Technology under NASA 
contract 1407. 
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Table 1. Flux densities of NGC 1377 in the Spitzer bands (this study ; see text) and in 

the ISO and IRAS bands (Roussel et al. 2003). The total photometric uncertainty for the 

Spitzer measurements is of the order of 10% in the IRAC bands, 5% at 24 /im, 10% at 

71 /im and 15% at 156 /xm ; for the ISOCAM measurements, ~ 20% ; and for the IRAS 

measurements, 10% at 12 and 25 /im and 15% at 60 and 100 /xm (IRAS Explanatory 

Supplement). 



band 


nuclear flux 


total flux 




(Jy) 




(Jy) 


3.6 /im (IRAC) 


0.0199^ 




0.0591 


4.5 /im (IRAC) 


0.0520 




0.0817 


5.7 /im (IRAC) 


0.2715 


= 


nuclear 


6.7/im (ISOCAM) 




0.366 


- 0.391*^ 


7.9 /im (IRAC) 


0.4105 


= 


nuclear 


12 /im (IRAS) 






0.43 


14.3 /im (ISOCAM) 




0.705 


- 0.721^ 


24 /im (MIPS) 






1.740 


25 /im (IRAS) 






1.83 


60 /im (IRAS) 






7.01 


71 /im (MIPS) 






5.64^ 


100 /xm (IRAS) 






5.97 


156 /im (MIPS) 






3.19^ 



'^If we try to decompose the surface brightness profile in the 3.6 /im band to estimate the 
central brightness of the disk, and adjust a point spread function to the residual central 
brightness above this level, we obtain an estimate of 0.0162 Jy for the nuclear flux, i.e. 10% 
lower. 

^The two flux values for each ISOCAM band correspond to measurements before and after 
correction for a slight saturation of the nuclear pixels. 

'^The MIPS 71 and 156 /^m fluxes include the latest flux calibration adjustments (+11% 
and +6% respectively, from version S12 to version S13 of the Spitzer Science Center pipeline). 
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Table 2. Line fluxes and upper limits and critical densities of the forbidden lines^ . All 

lines were measured at high spectral resolution, except H2 v=0 S(3) which is outside the 

wavelength coverage of the high-resolution modules. The uncertainties on the critical 

densities reflect a range of electronic temperature between 5 x 10^ K and 10"^ K. 



hne 


rest wavelength 


spectral order 


flux 


critical density 




(/xm) 




(10-i8Wm-2) 


(cm-3) 


H2 v=0 S(0) 


28.219 


LH14 


<9.15 




H2 v=0 S(l) 


17.035 


SH12 


65.43 ±2.94 




H2 v=0 S(2) 


12.279 


SH17 


39.14 ±10.03 




H2 v=0 S(3) 


9.665 


SLl 


65.10 ±7.64 




[Sill] 


34.814 


LHll 


< 11.93 


(1.5 ±0.3) X 10^ 


SIII] 


33.480 


LH12 


<7.20 


(6.1 ±0.9) X 10^ 


SIII] 


18.713 


SHll 


<2.00 


(2.0 ±0.3) X 10^ 


Nelll 


15.555 


SH13 


<8.97 


(2.4 ±0.4) X 10^ 


Nell 


12.814 


SH16 


17.38 ±2.13 


(6.0 ±1.0) X 10^ 



'^The critical densities are computed as: A^crit = -^ I \P (^c/lK)"°-^ ^(Tc) fi'u^]; 
where A is the transition probability, taken in the NIST atomic spectra database 
{http://physics.nist.gov/PhysRefData/ASD), C = 8.6287 x 10~^cm^s~\ and g^ is the sta- 
tistical weight of the upper level. We used collision strengths from Galavis et al. (1995) for 
the [SIII] hues and Dufton & Kingston (1991) for the [Sill] line. 
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(B— I) color with I— band contours . ^ 



-20° 55' 
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Fig. 1. — The (B-I) magnitude difference (with arbitrary zero-points because the maps are 
not ffux-cahbrated) is shown in grey scale, from to 1.2 mag. The maximum color in the 
center of NGC 1377 is 0.1 mag higher than displayed here. Isophotes of the I-band image 
are superposed for the intensities ui x 3", n = 1 to 5, where aj is the standard deviation of 
the sky brightness. 
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Fig. 2. — Residuals of tlie B-band image after subtraction of a smootli pliotometric model. To 
construct this model, radial brightness profiles were derived from small angular sections along 
the northern minor axis and along the eastern major axis, which appear undisturbed in the 
(B-I) color map. Then, the brightness distribution at any position angle 6 was interpolated 
by affecting a weight cos^^ (sin^^) to the major (minor) axis profile. Left: positive residuals. 
Right: negative residuals. Top: in units of the standard deviation of the sky brightness (in 
logarithmic scale for the negative component). Bottom: as a percentage of the local model 
intensities. The contours are as in Figure 1. 



3,6 jum 



4,5 jum 




Oa*^ 39™ 45^ 





1 ■ ■-'';,' ^v^"j^: ,= . 


. 




■^■ ° '". ■}, - -,.'■■ ' 




53 


&. ^ ' <^ ■ "•- ■*-. » 'JM) 


^ r * ^ - - 




.<^. ■-.' ■/ ^ -jjre 








»..^_. g^x ■ ■ 




^^^^^^^SC^J^ 


^^- ^ ^1; r , 


54 




^fe! 










'" "* * 1 ^ ^J ^ * 






--.^■•Of 


' - . 




. 4 L' 4fJL 


G- 


^ 5S- 


^ = .... VT : , 


t . " ^ "■' ^ - 




/c/ 


'"'. '' . * 






o' 




■'<}\i.H" i- 


p* 




' U- ''\ ^ • 1 






. * ' I" - 1 





03*^ 36™- 45^ 



03*^ 33™ 45^ 



Fig. 3. — Infrared images of NGC 1377 in the four IRAC bands and in the 24 /im MIPS 
band. The contours are the isophotes a x 3", n = 1 to 9 for all bands except 7.9 /im, and 
n = 1 to 10 for 7.9 /im. The structure of the images is explained in the text (Sect. 3). 
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Fig. 4. — Residuals of the 7.9 /xm image after subtraction of the mirrored northern half of 
the galaxy (see text). An intensity cut of 2 MJy/sr was applied. The brightest residuals have 
an intensity of ~ 100 MJy/sr, and the peak brightness of the nucleus is about 2000 MJy/sr. 
The contours are the B-band isophotes of intensities (Tb x (1 + 3"), n = 1 to 5. All the 
residuals are ghosts of the nucleus, except for the source labelled A. 
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Fig. 5. — Rest-frame optical spectrum in the 6500-6800A region, extracted in a 2.5" x 10" 
aperture centered on the nucleus. The vertical dotted lines indicate the wavelengths of the 
Ha, [Nil] 6583A, [SlI] 6716A and [SlI] 673lA hues. 
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Fig. 6. — Spectral energy distribution of NGC 1377 between (a) 3 and 40 /xm, in flux density 
units ; (b) 5 and 160 /xm, in power units and logarithmic scale. Notice the very deep 
absorption bands from amorphous silicates at 10 and 18 /im. Broadband imaging photometry 
from IRAC, ISOCAM (6.75 and 15 /im), MIPS and IRAS (60 and 100 /xm) is overplotted as 
grey diamond symbols with horizontal error bars representing the fllter widths. In (a), the 
spectra from different orders (SL2, SLl, LL2 and LLl in order of increasing wavelength) are 
separated by vertical dotted lines. In (b), the spectrum from a 90 K blackbody is overplotted 
as a grey dashed line. 
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Fig. 7. — H2 rotational lines in the mid- infrared spectrum of NGC 1377. Left: Rest-frame 
low-resolution spectra. Right: Rest-frame high-resolution spectra. The vertical dotted 
lines indicate the wavelengths of the transitions. The 12.29 //,m transition is detected in two 
different spectral orders at high resolution ; the two spectra are overplotted, shifted in flux 
for clarity. 
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Fig. 8. — Tentative PAH detections in the mid-infrared spectrum of NGC 1377. Left: Rest- 
frame low- resolution spectra. Right: Rest-frame high-resolution spectra. Different colors 
are used for different spectral orders. Spectra in the SH module were extracted in smaller 
apertures to increase the signal to noise ratio (in the central 2x2 pixels for the 11.25 /im 
feature and 4x4 pixels for the 17 /xm complex). For the 17 /xm complex, the vertical dotted 
lines indicate the wavelengths of the 17.1, 17.4 and 17.8 /xm features. The prominent 17.03 
H2 line is not labelled here and shown separately in Figure 7. In the middle and bottom 
spectra, the vertical lines indicate the wavelengths of the 11.3 /xm, 6.3 /xm and 6.7;um features. 
For each low-resolution spectrum, the spectrum of IRAS 08572+3915, normalized to that of 
NGC 1377 at 36 /im, is overplotted as a grey dashed line (see text). 
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Fig. 9. — Rest-frame high-resolution spectra of faint molecular absorption features and of 
the [Nell] line. Different colors are used for different spectral orders. The vertical dotted 
lines indicate the wavelength of the features. 
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Fig. 10. — The fit to the infrared spectral energy distribution of NGC 1377 by the DUSTY 
model, with the parameters discussed in the text, is overplotted in grey, for ry = 75, 85 and 
95. Broadband fiuxes are shown as big diamonds with their error bars. Color corrections 
between 1.00 and 1.08 have been applied to the far-infrared fiuxes (60, 71, 100 and 156 ;um) 
and corrections of 0.92 and 0.96 to the near-infrared fiuxes (3.6 and 4.5 /im). 
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Fig. 11. — Residuals after subtraction or normalization by the DUSTY model, (a) Minimum 
residual emission in the 5-15 /im range (for ry = 75). The peak wavelengths of the aromatic 
bands discussed in the text are shown by vertical dotted lines (at 6.3, 6.7, 7.7, 11.3 and 
12.7;um). We stress that due to modelling limitations, the 7.7 and 12.7/im bands are very 
uncertain (see text), (b) Residuals in the 15-35 /im range shown for ry = 75 (lower curve) 
and Ty = 95 (upper curve), in optical depth scale. The vertical lines indicate the positions of 
crystalhne sihcate features discussed by Spoon et al. (2006a) (at 16, 19.5, 23.5 and 28 /xm). 
The effect of fringing is seen in observations with the LLl module, particularly between 20.5 
and 24 /xm, and beyond 30 /xm. 
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Fig. 12. — Excitation diagram for the observed rotational and rovibrational transitions of 
H2. The abscissa is the temperature equivalent to the upper level energy of the transition, 
and the ordinate is the logarithm of the number of molecules in the upper level divided by 
the statistical weight, normalized by this ratio for the rotational S(l) line. The lower curve 
connects the data points uncorrected for extinction, and the upper curve corresponds to the 
maximum allowed extinction (see text), the normalization being done before the extinction 
correction. The v=l-0 S(0) upper limit, being inconsistent with the v=l-0 S(l) flux, was 
not considered, as it is probably underestimated. 
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Fig. 13. — Comparison of the rest-frame mid-infrared spectra of NGC 1377 (this study) and 
IRAS 08572+3915 (Spoon et ah 2006a). The spectrum of NGC 1377 is shown in black. The 
spectrum of IRAS 08572+3915, normalized to that of NGC 1377 at the long- wavelength end 
of the LLl module (36 /xm), is shown in grey. The far-infrared (40-120 /im) luminosity of 
IRAS 08572+3915 is ~ 5 x 10^^ Lg, about 100 times the far-infrared luminosity of NGC 1377, 
and it has a slightly higher -FeoAtrnZ-^iooMm ratio of 1.4. 
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Fig. 14. — Evolution of the luminosity in the 24 /xm MIPS filter as a function of redshift, for 
three different mid-infrared spectral energy distributions: those of NGC1377 (solid curve), 
M82 (dashed curve) and Arp220 (dot-dashed curve). M82 and Arp220 were chosen for this 
comparison because they are popular starburst templates. 



